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1 Introduction 

Weyl functions (also called Weyl-Titchmarsh or M-functions) and their gen- 
eralizations are an important and much used tool in the spectral theory 
of differential equations (see a necessarily small part of recent papers and 
books on this topic and various references therein [3-5, 10, 11, 14, 16-18,20, 
32,37,41,49,51]). Following the seminal work [27] and more general con- 
structions in [45,49] (see also some references therein), one can use struc- 
tured operators to solve inverse problems for Krein, Dirac, canonical, non- 
classical and non-self-adjoint systems and their discrete analogues (see, for 
instance, [2,15,16,37,40-42]). It was proved in [16,36,37,40-42] that the ker- 
nels of these structured operators are connected with the Weyl functions via 
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some kinds of Fourier transformations and can be recovered directly from the 
Weyl functions. It is also true that the Weyl functions can be represented as 
Fourier transformations of the corresponding kernels of the structured oper- 
ators. Important related works [18,50] on the A-amplitudes for Schrodinger 
operators gave rise to a whole series of interesting papers and results on the 
high energy asymptotics of Weyl functions and local Borg-Marchenko-type 
uniqueness theorems (see, for instance, [9,11,12,30] and references therein). 
See also further discussion of A- amplitudes in [5]. Finally, the interpola- 
tion of Weyl functions for scalar Schrodinger operators using the values of 
Weyl functions on a discrete countable set was dealt with in two interest- 
ing papers [6,35]. The second and more general paper [35] is based on the 
A-amplitude representation of Weyl functions. 

In this paper we consider self-adjoint Dirac systems and subclasses of 
canonical systems, which generalize Dirac systems. We consider direct and 
inverse problems, integral representations of Weyl functions, and interpola- 
tion of Weyl functions. Connections with Schur coefficients and A-amplitude 
are discussed. 

In Section 2 we formulate some previous results for Dirac systems, prove 
formula (2.21), which solves the inverse problem in a more general setting 
than previously considered, and derive for Dirac systems analogs of some A- 
amplitude-type results for Schrodinger equations. In Section 3 subclasses of 
canonical systems are studied and the case of rational Weyl matrix functions 
is discussed. Explicit solutions of direct and inverse problems are given for 
that case. 

A general case, where Weyl functions are not necessarily rational, is 
considered in Section 4. An integral representation of the Weyl function, a 
solution of the inverse problem for generalizations of Dirac systems directly 
via Weyl functions, and a local Borg-Marchenko-type theorem are obtained 
there. Section 5 is dedicated to interpolation of Weyl functions. 

As usual, by C we denote the complex plane, by C + we denote the open 
upper semi-plane, R is the real axis, R+ is the positive real axis, N is the 
set of positive integers, and No is the set of non-negative integers. We have 
(3' = -^(3. By Lp Xp (— oo, oo) we mean the space of pxp matrix functions with 
entries which belong to L 2 (— oo, oo). The pxp identity matrix is denoted 



by I p , the identity operator is denoted by I, and diag denotes a diagonal or 
block diagonal matrix. The set of bounded linear operators acting from Hi 
to H 2 is denoted by {Hi, H 2 } and a is used to denote the spectrum. 



2 Self-adjoint Dirac system: 

A-amplitude and inverse problem 

Consider the matrix self-adjoint Dirac system 

—u{x,z) = i{zj + jV(x))u(x, z), x > 0, 
where 
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(2.1) 



(2.2) 



I p is the p x p identity matrix, V is an m x m (m = Ip) matrix function 
with the p x p block entry v (often referred to as the potential) being locally 
summable. If v is summable on each finite interval [0, I], that is, iff is locally 
summable, there is (see Proposition 5.2 and its proof in [40]) a unique p x p 
matrix function ip(z) such that 



f°° ^ ^ r i 

/ I" /„ iip(z)* 1 u(x, z)*u(x, z) 

Jo L -^<f 



dx < oo (zeC: 



-icp(z) _ 
where u is the fundamental solution of (2.1) normalized by the condition 



(2.3) 



u(0,z) = K*, K:= — 



V2[I P ij 



v J 



(2.4) 



Definition 2.1 The p x p matrix junction tp(z) (z G C+), which satisfies 
(2.3), is called the Weyl function of system (2.1). 

The Weyl function of (2.1) is holomorphic in C + and Qtp > 0, which we 
discuss in detail in Section 5. 

We always assume that v is measurable. When v is also locally bounded, 
one can recover it from the spectral function or directly from the Weyl func- 
tion using Theorems 4.2 or 5.4 from [40], respectively. To recover v from the 



S(X) 



Weyl function one should notice that the Weyl function of system (2.1) on 
the semi-axis is a so called Weyl point and, according to Proposition 5.2 and 
formula (5.3) from [40], Theorem 5.4 from [40] is applicable. 

Let us describe the procedure suggested in Theorem 5.4 [40] to recover 
system on any intervals. For this purpose we need a family of operators 
which are associated with system (2.1) and have difference kernels: 

d f l 

Si = — — / s(x — t) ■ dt, s(x) = —s(—x)*, 
ax Jq 

where 

(2.5) 
Here (( + ir])~ 2 (p(( + iri) G L 2 xp (— oo, oo) for every fixed rj > 0, l.i.m. denotes 
the entrywise limit in the norm of L 2 (0, I) (/ > 0), and l.i.m. in the right- 
hand side of (2.5) is differentiable. Moreover, s(x) is boundedly differentiable, 
s(0) = |J P , and Si > for all / > 0. Thus, we have 

r l d 

S l = I + k(x-t) ■ dt>0, k(x) := s'(x) (s' = — s), k(x) = k(-x)*, 
Jo dx 

(2.6) 

and Si, Sf 1 E {^(0, 1), L 2 (0, 1)}. Notice, that for convenience purposes </?, s, 
and S in the present paper differ by scalar constant factors from the objects 
denoted by the same letters in [40]. 

We now introduce the matrix functions 

e 1 (x) = [i p o]u(x,o), e 2 (x) = [o i p ]u(x,o). (2.7) 

According to formula (4.16) from [40] we have 

9 2 (x) = ±= (\-I p I P ]-J X k(tyS^[2s(t) I p ]dty (2.8) 

where S' 1 (r = 2x) is applied to [2s(t) I p ] columnwise. It easily follows 
from (2.1) and (2.4) that 



u(x, 0)*ju(x, 0) = J, u(x, 0)Ju(x, 0)* = j, J :- 



I P 



I P 



(2.9) 



Using (2.1), (2.7), and the second relation in (2.9) we obtain 

v(x) = i6[(x)J6 2 (xy. (2.10) 

Finally, in view of (2.1), (2.4), (2.7), and the second relation in (2.9) we get 
equalities 

e 1 (o) = ^=[i p i p ], 9 1 (x)j9 2 (x)* = o, e' 1 (x)j9 1 (xy = o, (2.11) 

which uniquely determine 9\ assuming that 9 2 is already given. Rewrite 
Theorem 5.4 from [40] in a slightly modified form. 

Theorem 2.2 [40] Let tp be the Weyl function of the matrix Dirac system 
(2.1), where the potential v is locally bounded, that is, 

sup ||f(x)|| < oo for any / > 0. (2-12) 

0<x<l 

Thenv is uniquely recovered from tp by the formulas (2.10), (2.8), and (2.11), 
whereas s and S 2x in (2.8) are given by (2.5) and (2.6), respectively. 

Remark 2.3 Notice that according to formula (3.26) from [40] we have 

/•oo 

tp(z) = 2z e izx s(x)*dx, (2.13) 

Jo 

that is, the kernel s, which is used in Theorem 2.2 to solve the inverse prob- 
lem, is a Dirac system's analog of the A-amplitude. (To be more precise, 
k = s' is the analog.) Following M.G. Krein, the matrix functon k is called 
an accelerant in [2]. 

Clearly, formulas (2.5) and (2.13) are closely related (and (2.13) was used 
in [40] to derive (2.5).) Another closely related formula is representation 
(5.5) from [40]: 

POD 

<p(z) =z 2 e lzx X (x)dx, e- vx X (x) E L 2 vx J0, oo) for all r] > 0, (2.14) 



pxp\ 




where 



X (x) := -1% J s(t)*dt. (2.15) 

Jo 



According to p. 341 in [40] the operators Si admit a triangular factori- 
sation 

S t = E-^E- 1 )*, E l = I+ [ E(x,t) ■ dt, (2.16) 

Jo 

E' 1 = I + f F(x,t) ■ dt, sup (\\E(x,t)\\ + \\r{x,t)\\) <oo, (2.17) 

JO 0<t<x<l 

E~ 1 9 1 (x/2) = -^=[2s(x) Ip] (0<x<l), (2.18) 

v2 

where E h Ef 1 e {12(0,1), 12(0,1)} and Ef 1 is applied to 9\(x/2) colum- 
nwise. (Note that the expression for the triangular factor V of S^ 1 on p. 
341 in [40] is slightly more complicated but according to [44] we can put 
k(x) = l/v2 in that expression. Taking into account also that Si is here two 
times less than Si in [40] and so E[ is here \/2 times larger than its analog 
V in [40], we get precisely (2.16).) By (2.16) and (2.17) we have 

Sf 1 = E*E h I > 0, (2.19) 

and the matrix function E is locally bounded, that is, we can use equality 
(2.19) to determine pointwise the kernel of the integral operator in the left 
hand-side. Specifically, we get 

(Sf 1 f)(x)=f(x) + J (E(x,t) + J E(r,x)*E(r,t)dr)f(t)dt 

+ I (E(t,x)*+ f E(r,x)*E(r,t)dr)f(t)dt. (2.20) 

Now, the next result follows from Theorem 2.2 and formula (2.20). 

Theorem 2.4 Let ip be the Weyl function of the matrix Dirac system (2.1), 
where the potential v satisfies (2.12). Then v is uniquely recovered from ip 
by the formula 

v(L)=2t(S r 1 k)(l), (2.21) 

where k and Si are given by (2.5) and (2.6). 



Proof. By (2.19) we get S^ 1 = E* r E r . Note also that (E r f)(t) = (EJ)(t) 
for < t < r < I. Hence, formulas (2.8) and (2.18) imply 

e' 2 ( x ) = -V2(E l k\(2x)*(E l [2s I p ])(2x) = -2^,A;)(2x)*0i(x), 2x<l. 

(2.22) 

In view of (2.10) and the second equality in (2.11) one can see that 

0' 2 J0l = -0 2 J(0' 1 )* = -iv*. (2.23) 

By the first equality in (2.7) and the second equality in (2.9) we have 

X J0\ = I p . (2.24) 

From (2.22)-(2.24) we obtain v (x) = 2i(Eik)(2x) or, equivalently, 

v(^) =2i(Eik)(x), 0<x</. (2.25) 

In other words, for almost all values of x we have 

v(-) =2i(k{x)+ J E(x,t)k(t)dtY (2.26) 

Putting x = I in (2.20) and (2.26) we see that (2.21) is true. ■ 

Formulas related to (2.21) can be found in [2, 16] for the cases of continuous 
v and of a skew-self-adjoint Dirac system, respectively. 

The case of a continuous v was treated in [2] . If v is continuous, then by 
Theorem 1.2 in [2], the matrix function k is also continuous with a possible 
jump discontinuity at the origin. Moreover, it follows from the proof of 
Theorem 1.2 in [2] that the kernels E(x,t) and T(x,t) of the operators E\ 
and E^ , respectively, are continuous in the triangles < t < x < I. In 
particular, E(x,0) is well-defined. According to formulas (1.13) and (2.12) 
from [2] we have 

v(x) = -2iE(2x,0). (2.27) 

Moreover, since E(x,t) is continuous, formula (2.25) holds pointwise. By 
applying Ef 1 to both sides of (2.25) we get the following analog of the 
representation of A-amplitude of Schrodinger operator in [5]. (See formula 
(3.11) in [5], which is essential in the proof of Theorem 1 in the same paper.) 



Corollary 2.5 Let the potential v of system (2.1) be continuous. Then the 
matrix function k given by (2.5) and (2.6) admits the representation 



k(2x) = - % -(v{x) + 2 T(2x,2t)v(t)dt S 

= --(v{x)+2 I F(2x,2(x -t))v(x-t)dt). (2.28) 



The matrix functions k and T are continuous and (2.28) holds for all x G 

[0, oo). 

3 Canonical system. Explicit formulas 

3.1 Canonical system and Schur coefficients 

In this section we consider a canonical system 

—w(x,z)=izJH(x)w(x,z), H{x) = (3(x)*(3(x), x > 0, (3.1) 

dx 

where J is given by (2.9), /3 is an p x m (m = 2p) matrix function such that 

p(x) J/3(x)* = D = diag{di, d 2 , . . . , d p }, (3.2) 

and diag denotes a diagonal matrix. Weyl function of the canonical system 
on the semi-axis x>0 is apxp holomorphic matrix function tp(z), which 
satisfies the condition [49] 



[ I p iip(z)* ] w(x, z)*H(x)w(x, z) 



I P 



-i(p(z) 



dx < oo, z £ C 



+ ■. 



(3.3) 
where w is the fundamental solution of (3.1) normalized by w(0,z) = I2 P - 
The Dirac system of the previous section, is equivalent [40] to canonical 
system 

—w(x,z) = izJH(x)w(x,z), H(x) = u(x,0)*u(x, 0), x > 0. (3.4) 
dx 



Taking into account (2.9) one can see that a simple multiplication w(x, z) = 
e tzx w(x, z) of solution w by scalar factor transforms system (3.4) into system 
(3.1), where 



D = 2L 



v- 



p(x) = >/20i(z) 



(3.5) 



and 9\ is given by (2.7). Nevertheless, additional Weyl functions appear as 
a result of this transformation and the Weyl theory of system (3.1) requires 
some separate study even in this particular case. 

Before we turn to the general case (3.2), we consider the canonical system 
(3.1), (3.5) generated by Dirac system in greater detail. Separating j3 into 
two p x p blocks, namely j3i and (3 2 , and using (3.2) and (3.5) we obtain 



det f3 2 (x) ^ 0, ^{(3 2 {xT l (3i{x)) > 0. 
By (2.1), (2.7), and (2.9) we get 

f3'(x)Jf3(x)* = 26[(x)J6 1 (xy = 2iv(x)[0 I p }j 

According to (3.7) we have (3[f3^ + f3' 2 (3l = 0, that is, 

p' 1 {x) = -p l 2 {x)p l {xy(p 2 {xyy\ 



(3.6) 



(3.7) 



(3i 



It follows from (3.8) that {(3 2 l (5 l )' = - (5^ (5' 2 {(5 2 l (5 X + f3l((5* 2 )- 1 ) . Moreover, 
formulas (2.4), (2.7), and (3.5) imply f3 2 (0) = I p . Hence, we can recover (3 2 
from Z?^" 1 /?! using the following differential equation and initial condition 



fo = -/3 2 (/3 2 - 1 /3 1 )'(/3 2 - 1 / 3 1 + Pimr'Y 1 ' &(0) 



I p . 



(3.9) 



Remark 3.1 Note that fi 2 x fi\ is the continuous analog of the (well-known in 
the discrete case) Schur coefficients. Similar to the case of Toeplitz matrices, 
the function f3 2 1 f3 1 has the property ^R,(^f3 2 1 f3{j > 0. Moreover, system (3.1) 
corresponding to operators with difference kernels is uniquely recovered from 
fi 2 x fi\ via formulas (3.9) and fii = fi 2 ((3 2 l (3i) . 



3.2 Explicit formulas 



In the case of the rational Weyl functions and corresponding Hamiltonians H 
one can obtain explicit formulas for solutions of the direct and inverse prob- 
lems using a GBDT version of the Backlund-Darboux transformation [38,39]. 
See [8,13,19,22,24,31,38,52] and references therein for various versions of the 
Backlund-Darboux transformation and commutation methods. For explicit 
formulas via inversion of semiseparable operators see, for instance, [2,16,43]. 
Here we shall apply the procedure from [38] to the initial system 



d 



—w {x,z 

ax 



izJH w (x,z), Hq 



D/2 



L p j 



[ D/2 I p ] (det-D^O), 

(3.10) 
where D = D* = diag{di, d,2,---,d p } is a fixed diagonal matrix. (Notice 
that J in (3.1) and in (3.10) is slightly different from J in [38].) 

We consider a fixed integer n > 0, an n x n matrix a and two n x p 
matrices Ai(0) and A 2 (0), such that 

a-a* =iA(0)JA(0)*, A(0) := [Ai(0) A 2 (0)]. (3.11) 

We now introduce matrix functions A(x) and T,(x) with the equalities 
d 



dx 



A = -iaAJH , A(0) = [Ai(0) A 2 (0)]; 



(3.12) 



d 

dx 



S = *!**, E(0) = /„, * fc (x) := A(x) 



\k+i 



D/2 



k = l,2. 

(3.13) 

The first two equalities in (3.13) imply S(x) > /„. Thus, E is invertible and 
we can put 

H(x)=Vo(x)*H v (x), (3.14) 

where 

d 



-j-v = -q v , v (0)Jv (0y 



J, (3.15) 

q (x) := JA(x)*E(x)- 1 A(x)JH - JiJ ^A(x)*S(a;)- 1 A(a;). (3.16) 
The transfer matrix function, which we use in GBDT, has the form 

w a (x, z) = I m - iJA{x)*T l {x)- 1 {a - zI n )- 1 A{x). (3.17) 
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In view of the second equality in (3.10), the first equality in (3.12) and the 
third equality in (3.13), we see that (-^AjJA* = -ia*i**. Therefore, 
taking into account (3.11) and (3.13) we get 



— ( aS — S«* 
ax 



d 



■ i— UjA*) , «E(0) - S(0)a* = iA(0) JA(0)*. 

The two equalities above imply 

«E(x) - E(x)a* = iA(x)JA(x)*. (3.18) 

By direct calculation [45,49], it follows from (3.17) and (3.18) that 

w a (x,z)*Jw a (x,z) = J. (3.19) 

It follows from (3.1) and (3.15), respectively, that 

w(x, z)* Jw(x, z) = J, v(x)*Jv(x) — J. (3.20) 

According to Theorem 1 in [38] the fundamental solution w of the canonical 
system given by (3.1) and (3.14)-(3.16) is expressed via w a : 



w(x,z)—v (x) 1 w a (x,z)wo(x,z)w a (0,z) 1 t'o(0), 



(3.21) 



where Wq is the fundamental solution of the initial system (3.10). Here 
the fundamental solutions are normalized by the initial condition w(0,z) = 
w o (0,z) = I m (m = 2p). 

Remark 3.2 If det a^O, by formula (8) in [38] and by formula (3.15) we 
have 

v (x) = w a (x, 0)M, M := v (0), MJM* = J. (3.22) 

Moreover, for the case of the initial system (3.10) the matrix functions wq, 
A and S can be constructed explicitly. Hence, assuming det a ^ ; one gets 
explicit expressions for H (x) andw(x,z) via formulas (3.14), (3.17), (3.21), 
and (3.22). 
Indeed, put 



D/2 -D/2 



u 


'DO' 


, D = 




7 


_ 



(3.23) 
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It follows that 



Z- 1 = diagj/r 1 , D- 1 } 



D/2 
D/2 



(3.24) 



In view of (3.10), (3.23), (3.24) an d normalization condition wq(0,z) = I m , 
it is immediately apparent that 

w (x,z)=Ze izx °Z- 1 . (3.25) 

According to (3.12), to the second equality in (3.10) and to the third equality 
in (3.13), we have 

A^ = - ia y lD , ^(0) = A x (0) + ^A 2 (0) L>, * 2 (x) = Ai(0) - ^A 2 (0)L>. 

(3.26) 
The expression for \l/i follows directly from (3.26): 

\&i(x) = [exp(— idixa)fi exp(— id 2 xa)/ 2 ... exp(—id p xa)f p ], (3.27) 

[h f 2 ... f ,]:=*!$)■ (3-28) 

Formulas (3.26) and (3.27) provide explicit expressions for ^i and ty 2 - Fi- 
nally, from (3.13) and (3.23) we obtain 

A(x) = [* 1 (x) y 2 {x)}Z-\ X(x) = I n + f V 1 (t)V 1 (t)*dt. (3.29) 

Jo 

Definition 3.3 A canonical system given by (3.1) and (3.14), where H is 
defined in (3.10), vq is defined by (3.15) and (3.16), and A(x) and S(x) in 
(3.16) are defined by (3.12) and (3.13) or, eqmvalently, by (3.26)-(3.29), is 
said to be determined by the parameter matrices a and A(0) such that (3.11) 
holds. 

Further we always assume that (3.11) is valid. 

3.3 Direct problem: explicit solutions 

Remark 3.2 leads us to our next proposition. 
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Proposition 3.4 Let i>o(0) — I m (m — 2p). Then the matrix functions 



<p(z) = --D + * 1 (0)*(>y-zI r . 



1-1 



*v 



^) = -| J D|+§i(0)*(7-^/n)- 1 ^2, 



where 



7 =: a - i* 2 A 2 (0)*, ^i(O) = Ai(0) + -A 2 (0)L>, 

^2 = A!(0) - ^A 2 (0)£> = Ax(x) - ^A 2 (x)£>, 

7 =: a - i$ 2 A 2 (0)*, $x(0) = A x (0) - ^A 2 (0)|D|, 

§ 2 = Ai(0) + ^A 2 (0)|L>|, \D\ = diag{|di|, |d 2 |, . . . , |d p |} 



(3.30) 
(3.31) 

(3.32) 
(3.33) 
(3.34) 
(3.35) 



are Weyl functions of the canonical system determined by a and A(0). 

If D < 0, then the Weyl functions tp(z) and (p{z) coincide and, moreover, 
the Weyl function is unique. 

The matrices 7 and 7 given by (3.32) and (3.34), respectively, satisfy 
matrix identities 

7-7* = iA 2 (0)L>A 2 (0)*, 7* -7 = i($2-*i(0)) |^>| _1 (^ 2 — *i(0))*, (3.36) 

that is, for D < we obtain 17(7) l~l C + = and (p given by (3.30) does not 
have singularities in C+. 

Proof. The proof that tp is a Weyl function is somewhat similar to the proof 
of the corresponding fact for the self-adjoint Dirac system [21]. Put 



Q(z) :— w a (0, z)Z, Vt 2 {z) 



ui{z) 
w 2 (z) 



:= n(z) 







(3.37) 



where Uk are p x p blocks of Q 2 . In view of (3.21), (3.25) and (3.37) we have 



w(x, z)fl 2 (z) — Vq(x) 1 w a (x,z)Z 







L p J 



v (x) 1 w a (x,z) 



-D/2 ' 

(3.38) 
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Notice also that by (3.13) and (3.17) the equality 



-D/2 



Q{x) := [ D/2 I p ] w a (x,z 
holds. From (3.10), (3.14), (3.38), and (3.39) it follows that 



-i* 1 (x)*E(x)- 1 (a - zl n )- 1 * 2 

(3.39) 



Q 2 (z)*w(x, z)*H(x)w(x, z)tt 2 (z) = Q(x)*Q(x) (3.40) 

= %(a* - ^/ n )- 1 E(x)- 1 *i(x)*i(x)*E(x)- 1 (a - zl n y l ^ 2 . 



According to (3.13) we have 
d 



dx 



E- 1 = -E-^x^E- 1 , E(0) = I n , E(x) > I n - 



(3.41) 



Hence, from (3.40) it follows that 

Q 2 (z)*w{x, z)*H(x)w(x, z)Vt 2 (z)dx = %(a* - zQ -1 (3.42) 

'o 

x (E(O)- 1 - E(r)- 1 ) (a - zI n )-^ 2 < %(a* - zl n )~\a - zI n )-^ 2 . 
Compare (3.3) and (3.42) to see that the function 

ip(z) = iu 2 (z)ui(z)~ l , (3.43) 

where Uk are the blocks of Q 2 , satisfies (3.3) (excluding, possibly, a finite 
number of points), i.e., tp of the form (3.43) is a Weyl function of system 
(3.1), (3.14). 

Next, let us show that the right-hand sides of (3.30) and (3.43) coincide. 
Since i>o(0) = I m and E(0) = I n , we use (3.17), (3.23), and the third equality 
in (3.13) to rewrite (3.37) in the form 

u^z) = I P - tA 2 (0)*(a - zl n )'^ 2 , u 2 (z) = ~D - iAx(0)*(a - zI n )-^ 2 . 

(3.44) 
The following procedure is a standard one in system theory. Rewrite the first 
equality in (3.32) as 



i^ 2 A 2 (0)* = a - 7 = (a - zl p ) - (7 - zl p ) 



(3.45) 
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to check that 

(I p - iA 2 (0)*(a - zl n )- 1 ^ 2 )(I p + *A 2 (0)*( 7 - zl n )' 1 ^ 2 ) = I p . 
In other words, we have 

co^z)- 1 = I P + 2A 2 (0)*( 7 - zI n )'H 2 . (3.46) 

From (3.44)-(3.46) it follows that 

(3.47) 



1<j0 2 \Z)IjO\\Z) 



--D + ^ 1 (oy( 1 -zi P r 1 ^ 2 , 



that is, the function given by (3.30) is a Weyl function. 

Now, consider the matrix function ip, given by (3.31). If D < 0, we have 
tp — ip and the case of tp is as above. If the inequality D < does not hold, 
we can introduce two diagonal matrices Pj = P* = P 2 (j — 1, 2) such that 



P 1 + P 2 = I P , D(P 1 -P 2 ) = \D\. 



(3.48) 



Here Pi has nonzero entry (entry equal to 1), when the entry of D in the 
same row and column is positive, and P 2 has nonzero entry (entry equal to 
1), when the corresponding entry of D is negative. Let us show that for 
Qz > \\a\\ we have 



[Pi P 2 ]n(z)*w(x, z)*H(x)w(x, z)Q(z) 
It suffices to show that 

[0 P 2 ]tt(zyw(x,z)*H(x)w(x,z)n(z) 



Pi 
p 2 





p 2 



dx < oo. 



dx < oo, 



(3.49) 



(3.50) 



which is immediately evident from (3.42), and that 



[Pi ti\tt(zyw{x,zyH(x)w{x,z)tt{z) 



Pi 





dx < oo. 



(3.51) 



To prove (3.51) note that similar to (3.40) one obtains 



[Pi ti\n(zy w {x,zyH(x)w{x,z)tt{z) 



Pi 





Qi(x)*Qi(x), (3.52) 
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Q x {x) = (D- i*i(x)*E(x) _1 (a - zI n )- 1 ^f 1 (x))e izxD P 1 . (3.53) 

It is an immediate consequence of the above that the entries of De lzxD Pi 
belong L 2 (0, oo). It follows from (3.27) that the entries of ^ \{x)e lzxD P\ 
are bounded for Qz > \\a\\. Because of (3.41) the entries of 1 ifi(x)*'E(x)~ 1 
belong L 2 (0, oo). Thus, in view of (3.53) the entries of Qi belong L 2 (0, oo) 
(^sz > || a ||), and using (3.52) we see that (3.51) is true. Hence (3.49) is also 
valid. Substitute (p of the form 



<p(z) 



IU 2 {Z)UJ 1 {Z 



,-1 



Ui{z) 

v 2 (z) 



Q(z) 



Pi 
Pi 



(3.54) 



into the left-hand side of (3.3). By (3.49) the function (p of the form (3.54) 
satisfies (3.3) in the domain Qz > \\a\\. 

According to (3.37), (3.48), and (3.54) we have 



bJl(z) 

£2(2) 
In view of (3.55), the formula 



w a (0,z) 



\D\/2 



(3.55) 



lU 2 [Z)UJi{Z) 



|D| + ^i(0)*(t - ^/„) _1 *2 



(3.56) 



can be proven quite similar to (3.47). Equation (3.31) is an immediate con- 
sequence of (3.54) and (3.56), and so <p given by (3.31) satisfies (3.3). 
To prove the uniqueness of the Weyl function for D < 0, notice that 



Hq + cJ 



D/2 
2D- l {\D + cI p ) 



AcD- 1 {D + cI p )D- 1 



[D/2 2{\D + cI p )D- 1 ] 



(3.57) 



When D < and the scalar c > is sufficiently small, formula (3.57) implies 
H > —cJ. Hence, because of (3.14)-(3.16) we have 



H(x) > —cv (x)*Jv (x) = —cJ. 



(3.5* 
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Using (3.1) we obtain 

— c— (w(x, z)*Jw(x, z)) = ic(z — z)w(x, z)*H(x)w(x, z) > 0, z E C + . 

ax v y 

(3.59) 

From (3.58) and (3.59) it follows that 

w(x, z)*H(x)w(x, z) > —cw(x,z)*Jw(x,z) > —cw(0,z)*Jw(0,z) = —cJ. 

(3.60) 
Formula (3.60) implies that 

r r i i 

/ [ I p —I p ]w(x,z)*H(x)w(x,z) p dx>2crl p . (3.61) 

Jo ' L — *p J 

Denote by L C C m (m = 2p) the subspace of vectors / such that 

f*w(x,z)*H(x)w(x,z)fdx < oo. 



Inequality (3.61) implies dimL < p, which implies the uniqueness of the 
Weyl function. 

Finally formulas (3.11) and (3.32)-(3.35) imply (3.36). ■ 

Remark 3.5 Representations like (3.30) and (3.32) are called realizations 
in Control Theory. 

Remark 3.6 It follows from the proof of Proposition 3.4 that ip satisfies 
(3.3) in C + \ (o~(a) U (7(7)) and (p satisfies (3.3) in {z : Qz > ||a||}. 

3.4 Inverse problem: explicit solution 

The following theorem allows us to recover H explicitly from the rational 
Weyl functions. 

Theorem 3.7 Let the diagonal p x p matrix D be negative (i.e., let D < 0). 
If tp is a rational matrix function such that 

Km (p(z) = Ud\, %ip(z)>0 (zeC+), (3.62) 
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then if admits a realization 

<p(z) = l -\D\ + *i(0)*( 7 - zl n )~ 1 y 2 , (3.63) 

where ^i(O) and \I/ 2 are n x p matricies, nGN, and nx n matrix 7 satisfies 
the matrix identity 

7 - 7* = 1(^(0) - ^D" 1 ^!^) - * 2 )*. (3.64) 

Moreover (p is the Weyl function of the canonical system determined (in the 
sense of Definition 3. 3) by the parameter matrices 

a = 7 + ^ 2 A 2 (0)*, A 1 (0) = ^(V 1 (0) + V 2 ), A 2 (0) = (q 1 (0)-q 2 )D- 1 . 

(3.65) 

Proof. Realization (3.63) follows directly from Proposition 4.1 in [43] (see 
also Theorem 5.2 in [22]). The identity (3.11) follows from (3.64) and (3.65), 
that is, the requirement (3.11) for parameter matrices is fulfilled. Moreover, 
(3.65) implies relations (3.32) and (3.33). Now, compare (3.30) and (3.63) to 
see that tp is the Weyl function by Proposition 3.4. ■ 



4 Canonical system. General formulas 

In the Krein paper [27] (see also [1,2]) Krein system was treated as a system 
generated by an accelerant k(x) = k(—x)*, where k was the kernel of operator 
S (with difference kernel) of the form (2.6). In this section we assume that 
the diagonal matrix D is negative (i.e., D < 0). In a similar way to [27], 
operators with \D [-difference kernel generate other subclasses of canonical 
systems (see [43,48,49] and references therein): 

Proposition 4.1 Let k(x) = {kij(x)}?j =1 be a p x p matrix function such 
that 

k(x) e L 2 pxp (0,l), k(-x) = k( X y, S h S^e{Ll(0,l),Ll(0,l)}, (4.1) 

S t :=I+ / {kij(djt - dix)} p ij=1 ■ dt > 0, for all < / < 00, (4.2) 



that is, operators Si determined via k are positive, bounded, and boundedly 
invertible. Then operators Si admit a triangular factorization 

S^ 1 = E*E h Ei = I+ [ E(x,t) ■ dt G {Lj(0, I), L 2 (0, I)}, (4.3) 

Jo 

where E(x,t) is a Hilbert- Schmidt kernel. 

Moreover, k generates, in terms of S and E, the canonical system (3.1) 
such that (3.2) holds, and the Hamiltonian H = (3*f3 of the canonical system 
is given via 

f3(x) = (E t U)(x), U(x) := [D{s l3 {\d,\x)}l J=1 J p ] , (4.4) 

1 f x 

s(x):=-I p +\D\- 1 k(t)dt. (4.5) 

2 Jo 

This H(x) is summable on all finite intervals (0, I). 

The fundamental solution of the canonical system is given by 

w(l, z) = I 2p + izmiSr\l - zA)- l Yi h A, Hi e {C 2p , Lj(0, /)}; (4.6) 

A = Ai = iD f -dt; Uig = U(x)g, g G C 2p , (4.7) 

Jo 

where C 2p is the 2p- dimensional vector space. 

Proof. According to [46] (Ch. 6) the operators A, Si, Hi given by (4.2) and 
(4.7) satisfy the operator identity 

AS t - S L A* = iliiJIil (4.8) 

that is, they form an S'-node [45,46,48,49]. By Theorem 4.2.1 [49] there is 
a system corresponding to the family A, Si, Hi (0 < I < oo) of S-nodes, and 
w given by (4.6) is its fundamental solution. Moreover, if B{1) := n*^ -1 !^ 
is differentiable, this system is canonical, and its Hamiltonian is given by 

"<" = ><" = M n ' s - ln - (49) 

The factorisation ST" = (I + E + )(I + E_) , where E + (EJ) is an upper (lower) 
triangular integral operator, is clear from the factorization result on p. 184 
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in [23]. Now, taking into account Si = S* it is easy to derive (4.3) (see, for 
instance, Section 4 in [43]). According to (4.3) and (4.9) we get H = (3* (3, 
where /3 is given by (4.4). Note that E(x,t) does not depend on I since the 
lower-upper factorisation of Si is unique, and the lower-upper factorisation 
of Si is unique because the lower-upper factorisation of / is unique. 

It remains to be shown that the equality (3.2) holds for almost all values 
of x. For that purpose we multiply (4.8) by E\ from the left and by E\ from 
the right. Then, taking into account (4.3) and (4.4) we get: 

EiAEf 1 - (EiAEf 1 )* = i(EiYi)J(EiY[)* = if3(x)J f (5(t)* ■ dt. (4.10) 

Jo 

As EiAEf 1 is a lower triangular operator, formula (4.10) implies 

EiAEf 1 = i(3(x)J [ p(t)* ■ dt. (4.11) 



Denote the kernel of E l 1 by E and rewrite the equality (4.11) in terms of 
the kernels of integral operators: 

/x rx ry 

(DE(y,t) + E(x,y)D)dy+ E(x,y)D E(r,t)drdy (4.12) 

= /3(x)J/3(t)* (t < x). 

Equality (3.2) follows from (4.12). ■ 

Definition 4.2 The class of Hamiltonians H(x) (0 < x < oo) of canoni- 
cal systems generated by matrix functions k, which satisfy the conditions of 
Proposition 4-1, is denoted by %{D). 

We next consider some fixed H € %{D) and turn our focus towards Mobius 
(also called linear-fractional) transformations 

y{z, I) = t(W 11 (z)V 1 (z) + W 12 (z)V 2 (z)) {W 21 (z)V 1 (z) + W 22 (z)V 2 (z)y\ 

(4.13) 

where 

W(l,z) = {W -(*)}?, =1 = w(l,zY, (4.14) 
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and the pairs {Vi(z), Vi(z)) are pairs ofpxp matrix functions, which are 
meromorphic in C + and satisfy (excluding, possibly, a discrete set of points) 
the following relations 

v\Vx + v* 2 v 2 > o, v{v 2 + v;v x > o, (4.15) 

that is, {Vi, V2} are pairs of nonsingular matrix functions with property- 
J. The set of Mobius transformations (4.13), where the pairs {Vi, V2} vary 
and I is fixed, is denoted by A/]. It follows from Statement 3 in [36] and 
interpolation results in [45] (see also their formulation in Theorem 1 from [36] 
or Chapter 1 in [47]) that there is a unique matrix function 

<p=c\m, (4.16) 

l>0 

and this tp admits representation 

/>oo 

<p(z) = zD 2 l e-^ D {^(|d,|x)}f J=1 dx (ze<C+), (4.17) 

where s is given by (4.5). Moreover, according to Theorem 1 from [36] the 
scalar products {Sif, f) (0 < / < 00) admit representation 

(y e- UxD f(x)d x ydr(t)( y j e-^ D f(y)dy). (4.18) 

Here r is the p x p non-decreasing matrix function from Herglotz represen- 
tation of ip: 



f°° ( 1 t \ f c 



l + tz dr(t) 
t-z 1 + t 2 ' 
(4.19) 
where p, > (p — in our special case), 1/ = v*, and 

/OO 
(l + t 2 )- 1 dr(t)<cx). (4.20) 

•OO 

Note that according to (3.1) and (4.14) we have 
— (w(x,z)*Jw(x,z))=0, andso W(l, z) = Jw(l, z)' 1 J. (4.21) 
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From (4.21) we get 

(W(l,z)*)~ 1 JW(l,z)- 1 = Jw(l,z)*Jw(l,z)J. (4.22) 

Moreover, taking into account (4.6) and (4.8) we obtain 

w{l, z)*Jw(l, z) = J + i{z- z)Ii*(I - zA^Sf 1 ^ - zA)- 1 ^. (4.23) 

Using (4.4) and (4.7), separate 11; into two blocks $ fe G {C p , Lj(0, /)}: 

$ih = D{ Sij (\di\x)} p id=1 h, ^ 2 h = h (heC p ); n, = [$i $ 2 ]. (4.24) 

We omit index "/" in A and $& as the choice of "Z" is clear from the context. 
Finally, formulas (4.22)-(4.24) imply 



r(l, z) : = -[I p 0] (W(l, z)*) 1 JW(l, z)' 1 



h 



(4.25) 



i(z - z)®* 2 (I - ■zA*)- 1 S~ 1 (I - zA)- 1 ^ > (ze 



By (4.16) and (4.25) the conditions of Proposition 9.1.5 from [49] are fulfilled, 
and hence tp satisfies (3.3). 

Theorem 4.3 Let H G %{D) (D < 0). Then the unique Weyl function of 
the corresponding canonical system (3.1) admits the representation 

/•oo 

<p(z) = -zD / e izx s(x)dx (z G C+), (4.26) 

Jo 

where s is given by (4.5). 

Proof. Representation (4.26) is clear from (4.17). 

It remains only to prove that ip satisfying (3.3) is unique. Here, it would 
be of interest to give a proof connected with S'-nodes to compare with the 
proof of uniqueness from Proposition 3.4. For this purpose let us consider 
$*(/ - zA^Sf^I - zA)- 1 ^. 

First, fix some a G 1R + , put I = na (n G N), and consider the scalar 
product in (4.18). It is clear that 

v 
W/l^EW^'M), /=:{/;}?=!, e, :={^}f =1 . (4.27) 
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Formula (4.18) and the equality below 

rna ra pa 

Jo k=1 Jo 

/^-/.(^(fc-ljajeL^a), 

imply that 

n /'oo pa 2 

{SnafjGjJjej) <nJ2 / e ~ Uydj fjk(y)dy e*dr{t)ej. 

k=l J -°° J° 

It follows from (4.18) and (4.29) that 

n 
(SnafjGjJjej) < n\\S a \\ ^ WfjkW 2 = Tl\\ S a \\ \\ fj \\ 2 . 



' jk (y)dyej, (4.28) 



(4.29) 



(4.30) 



fc=i 



and so we get 
-i 



By (4.27) and (4.30) we have (S na f,f)<pn\\S a \\\\f\\ 2 , 

\\Sna\\ <pn\\S a \\, i.e., S~l > (pn\\S a \\) ± I. (4.31) 

One can easily check that (/ — zA)~ 1 § 2 — e lzxD . Thus, using (4.31) we derive 

mi-zA^S^il -zA)- 1 <S> 2 >c(n)I p , c(n) -> +oo (n -> oo). (4.32) 

Now, let ip and ip, such that ip(z) := ip(z) — (fi(z) ^ for some z G C+, 
satisfy (3.3). This implies that 



[ iip(z)* ] w(x,z)*H(x)w(x,z) 

According to (3.1) and similar to (3.59) we have 
d 





-iijj(z) 



dx < oo. 



(4.33) 



dx 



w(x, z)*Jw(x, z) J — i(z — z)w(x, z)*H(x)w(x, z). 



Hence, we rewrite (4.33) as 



sup [i(z — z) 1 [0 iip(z)* ]w(l,z)*Jw(l,z) 

l<OCl 





—iip(z) 



< oo. (4.34) 
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In view of (4.23) inequality (4.34) implies 

sup U(z)*mi - zAT'S-^I - zAy^iz)) < oo, 

which contradicts (4.32). ■ 

Theorem 4.3 yields a solution of the inverse problem directly via Weyl func- 
tion. 

Theorem 4.4 Let (p be the Weyl function of the canonical system such that 
its Hamiltonian H G 'H(D) (D < 0). Then we have 

e-^s(x) G (Lj xp (0, oo) f) L 2 pxp (0, oo)) (4.35) 

for any i] > 0, and the matrix function k(x) (x > 0) is recovered via Fourier 
transform 



k(x) 



^. u ,,_/y.^ 4 (4 , 6) 



The Hamiltonian H is recovered from k as in Proposition J^.l or, equivalently, 
using (4.9). 

Proof. To prove (4.35) recall the first inequality in (4.31), which holds for 
all a G K + . It follows that 

IN| < ||S- rri || <p(l + 1)115x11 <Ci(J + l), Ci ER+. (4.37) 

We next represent the operator $i = $u given by (4.24) in the form 

$x = T l {z) + i$ 2 <p(z) = (J - zA)(I - zAy'T^z) + t<S> 2V (z), (4.38) 

Now, we need inequality (22) from [36]: 

||(7 - zAy'T^W < \\S l \\^^/\\(<p(z)-<p(zy)/(z-z)\l ^ GC+ . (4.39) 

Note that we can fix in (4.38) and (4.39) any z G C+. Taking into account 
relations ||$2|| = Vt and ||A|| < Cy, we derive from (4.37)-(4.39) that 

||$i|| <C 3 (/ + 1) 3/2 . (4.40) 
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Relation (4.35) follows from (4.24) and (4.40). 
Rewrite (4.26) as 

poo 

\D\- 1 <p({ + irj)/({ + iri)= e Kx e'^ s{x)dx (C e R, rj > 0). (4.41) 

Jo 

In view of the Plancherel theorem formulas (4.35) and (4.41) imply 

|£>|-V(C + iv)/(C + iv) = l-i-m.a-.oo / e iCx e^ x s{x)dx, (4.42) 

Jo 

e-^six) = ^-l.i.m. _ HX) f e-^l^rVlC + m?)/(C + mjK (4-43) 

for all fixed 77 > 0. Note that (4.43) holds for x > 0, that is, l.i.m. in (4.43) 
is considered in Lp Xp (0, 00). Finally, equalities (4.5) and (4.43) yield (4.36). 
■ 

Using Theorems 4.3 and 4.4 we obtain our next Borg-Marchenko-type result. 

Theorem 4.5 Let tp and <p be the Weyl functions of the canonical systems 
with Hamiltonians H and H , respectively, where H, H G %{D) (D < 0). 
Let the equalities 

e- iza (ip(z) - <p(z)) = 0(1) (z -»■ 00) (4.44) 

hold for all a < I on some ray Qz/tflz = c in C+. Then we have 

H(x) = H(x) for 0<x<l/d, d := max \d k \. (4.45) 

Proof. We mark functions corresponding to tp with a tilde, e.g., H, 's. 
Consider the entire matrix function 



a 



u{z):= e lz ^~ a \s(x)-^(x))dx. (4.46) 

Jo 

It is clear that ||o;(z)|| is bounded in the closed lower semi-plane and tends 
to zero on some ray Qz/^tz = C\ there. To estimate u on the ray ^sz/tflz = c 
in C + we use (4.26) and write 



u(z) = e~ iza {z\D\)- 1 (cp(z) - tp{z)) - / e iz( - x ~ a) (s(x) - s{x))dx. (4.47) 



DO 



II 
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From (4.35), (4.44), and (4.47) one gets the boundedness of ||a;(z)|| on the 
ray ^sz/$tz = c. Hence, according to the Phragmen-Lindelof theorem ||a;(z)|| 
is bounded in C + , and thus also in C. As ||o;(2;)|| tends also to zero on some 
ray, we see that oo(z) = 0, that is, s(x) = li(x) for x < a. Since this fact is 
true for all a < I the equalities 

s(x) = s(x), k(x) = k(x) (0 < x < I) (4.48) 

follow. Now, take into account (4.2), (4.4), (4.9), and (4.48) to obtain (4.45). 



5 Interpolation of the Weyl function 

As can also be seen from the cases treated in Sections 2-4 Weyl functions 
contain important information about systems, and their interpolation is of 
particular interest. One of the possible approaches is interpolation by rational 
Weyl functions. Recall that an inverse problem for rational Weyl functions 
of a subclass of canonical systems was solved explicitly in Subsection 3.4. In 
this section we shall consider another approach, namely an approach in the 
spirit of [35]. Theorem 4 from [35] (after reformulation for the matrix case 
and for the upper semi-plane) has the following form. 

Theorem 5.1 Let matrix function F admit the representation 

poo 

F(z) = / e izx f(x)dx, e'^fix) E Lj xp (0, oo) for all r] > 0. (5.49) 
Jo 

Then for any e > and Qz > | + e we have 

oo . n 

F(z) = ^2c n (z+ - -ie)^2a nq F(iq + ie), (5.50) 

n=0 q=0 



where 



mn + q)l c w (A)^ (2 ^ )n ^ ( fi^ A) . (5.51) 



(q\y(n-q)\' ny ' n^(?+i-iA) 
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The series in (5.50) converges uniformly on any compact subset in the semi- 
plane Qz > | + e and we get 

N . n 

\\F{z) -Y J c n {z+ l --ie)Y j a nq F{iq + ie)\\ = o(ivl +£ - 32 ), TV -+ oo. 

n=0 <?=0 

(5.52) 

Our next theorem is an analogue for the Dirac system case of the main 
Theorem 5 from [35]. 

Theorem 5.2 Let tp be the Weyl function of system (2.1), where v satisfies 
(2.12). Then for any e > and Qz > \ + £ the matrix function tp admits the 
representation 

oo . n 

tp(z) = -z 2 ^c n (z + - - is)^2(q + ey 2 a nq >p(iq + ie). (5.53) 

n=0 <j=0 

Moreover, for Qz > | + s we get 

N . n 

\\v( z ) + z2 ^2 c n{z + -~ is) ^(q + ey 2 a nq ip(iq + ie)\\ 

n=0 q=0 

= 0(n^ +£ -* z ), N^oo. (5.54) 



Proof. The theorem's statement follows from formula (2.14) and Theorem 
5.1. ■ 

Remark 5.3 One can use various transformations of tp to change the inter- 
polation set {z : z = iq + ie, q G N }. The simplest example is tp(z + z ), 
where z G (C+ U R). Noticing that tp(£)jz 2 satisfies the conditions of The- 
orem 5.1, we see that tp(z + Zq)/(z + z ) 2 (as a function of z) also satisfies 
these conditions. Hence, we obtain 



tp(z + zq) = - (z + z ) 2 y^ y c n (z+ - - ie) 



x ^2(q + £ - iz ) 2 a nq ip(z + iq + ie). (5.55) 

<?=o 
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Clearly, representation (4.26) enables us to apply Theorem 5.1 to Weyl 
functions of canonical systems. 

In fact, interpolation formulas of (5.53) type are true for the Weyl func- 
tions of the much wider class of Dirac systems as well as for other classical 
and non-classical systems because the corresponding Weyl functions are ei- 
ther bounded or belong to the Nevanlinna (also called Herglotz) class, that 
is, Qip(z) > 0. 

Indeed, according to Proposition 4.2 in [40] there is a unique Weyl func- 
tion f associated with a Dirac system on [0, oo) with a locally summable 
potential v. Moreover, this f> belongs to the Nevanlinna class and therefore 
admits a Herglotz representation (4.19). 

Proposition 5.4 Suppose (p(z) (z G C+) belongs to the Nevanlinna class. 
Then the matrix functions [z + i5)~ 2 (p(z + iS) (5 > 0) satisfy the conditions 
of Theorem 5.1. 

Proof. Put z + i5 = £ + irj (£ G R, rj > 5). Then, in view of (4.19) and 
(4.20) (and for some C\, C<i > 0) we get the inequality 
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\^ + irj\- 4 \h*<p^ + iv)h\ 2 dq 

oo / 
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where /i G C p , /i*/i < 1, T/i(t) = h*r(t)h, and C2 does not depend on h and on 
rj > 5. Therefore, the norms of entries of (£ + ii])~ 2 Lp(l; + ^77) in L 2 (— 00, 00) 
are uniformly bounded for rj > 5. Thus, we can apply Theorem V from [33] 
and derive 

(z + i5y 2 ip(z + iS) = Li.rn.a_oo / e izx f(x)dx ($$z > 0), (5.57) 

Jo 
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where / G Lp Xp (0, oo). Representation (5.49) follows from (5.57) and from 
the fact that (z + i5)~ 2 (p(z + iS) is analytic. ■ 

Remark 5.5 The integral representation (5.49) is true also for functions tp 
from the generalized Nevanlinna class N K under additional conditions 

lim tp(ir)) = 0, lim rj\^>ip{ir])\ < oo 

in some semi-planes Qz > h v [28]. Weyl functions from the generalized 
Nevanlinna class appear, for instance, in [25, 26, 29, 34] (see also references 
therein). 
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